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We present measurements of the lifetimes and polarization amplitudes for B® — > J/ip4> an d 
— ► J/ip K*° decays. Lifetimes of the heavy (H) and light (L) mass eigenstates in the B® 
system are separately measured for the first time by determining the relative contributions of 
amplitudes with definite CP as a function of the decay time. Using 203 ± 15 B® decays we 
obtain t l = (1.05 to ll ± 0.02) ps and t h = (2.07 ±^46 ± °- 03 ) P s - Expressed in terms of 
the difference Ar s and average P s , of the decay rates of the two eigenstates, the results are 
Ar s /r s = (65 I 2 ! ± 1)%, and Ar s = (0.47 to ll ± 0.01) ps" 1 . 

DOI: 10. 1103/PhysRevLett. 94. 101803 PACS numbers: 13.25.Hw,11.30.Er,14.40.Nd 



Particle-antiparticle oscillation occurs for both and 
B® mesons and gives rise, in each system, to two eigen- 
states with definite masses (heavy, to# and light, m^) 
and widths (Th and Tl)- In the Standard Model 



(SM) , this oscillation is due to second-order contributions 
from the weak interaction and depends on the Cabibbo- 
Kobayashi-Maskawa (CKM) quark mixing matrix. Os- 
cillation has been observed in the B® system, and the 



4 



mass difference (Am = m# — mj,) is Am^ = (0.507 ± 
0.007) ps _1 l]. In the B® system, direct observation of 
the oscillation signal has been a challenge: the 95% CL 
limit for the mass difference is Am s > 14.4 ps _1 The 
ratio of the decay width difference Ar = Tl — T# to the 
average decay width T = (Tl + r#)/2 is expected to be 
small (0.2% - 0.3%) for the B° d system @, but sizable for 
the B° s system 0. Ref. predicts AT S /T S = (12 ±6)% 
and Ref. [j| gives the ratio of the decay width difference 
to the mass difference. If Am s is too large to be directly 
measured, a measurement of AT S could serve instead, 
along with Amj, in tests of the unitarity of the CKM 
matrix. In the SM, the mass eigenstates in the B° s sys- 
tem are expected to be nearly CP eigenstates. The light 
mass eigenstate is expected to be CP-even and to have 
a larger decay width, and thus a shorter lifetime, than 
the heavy mass eigenstate . By exploiting decays with 
known CP content, it is possible to measure the two de- 
cay widths separately. 

The decays B° -> J/tP4> and B^ -> J/tpK*°(892) are 
pseudoscalar to vector-vector transitions and are char- 
acterized by three amplitudes. These amplitudes corre- 
spond to transitions in which the J/tp and <j> (or K*°) 
have a relative orbital angular momentum L of 0, 1, or 
2. In the transversity basis 0, the decay amplitudes 
correspond to linear polarization states of the vector 
mesons. The L = 1 decays take place via the decay 
amplitude A± and correspond to a parity-odd (perpen- 
dicular) correlation between the transverse linear polar- 
ization states of the vector mesons. The other two decay 
amplitudes Aq and Au lead to decays corresponding to 
linear combinations of the parity-even L = and L = 2 
decays. In this analysis, the fully reconstructed decays 
B° s -> J/tpcp (with J/ip -> fj+fi- and 4> -> K+K-) and 
B a d -> J/iPK* (with J/tp -► n + n~ and K*° -> K + -k~) 
and their charge conjugates are used to measure the po- 
larization amplitudes. The observed final state particles 
for the B° and B s decays K + K~) have a definite 

CP, which depends on L, and a definite angular distri- 
bution. We determine the decay widths for the heavy 
and light B® mass eigenstates by measuring the relative 
contribution of the CP-odd and CP-even decays to the 
observed angular distribution as a function of the decay 
time. The B d decays provide a valuable control sample 
since they are expected to occur via similar (parity-odd 
and parity-even) decay amplitudes Q. 

The analysis uses a portion of the data from Run II 
at the Fermilab Tevatron pp collider, corresponding to 
an integrated luminosity of about 260 pb -1 . The data 
were collected with the upgraded Collider Detector at 
Fermilab (CDF) 0, the most relevant components of 
which are described below. A five-layer double-sided sil- 
icon microstrip detector, SVX, provides track measure- 
ments at radii between 2.5 and 10.6 cm and allows for pre- 
cise vertex reconstruction. A cylindrical drift chamber, 
COT, with eight alternating axial and stereo superlay- 



ers (each super-layer containing 12 sense wires) provides 
track measurements for charged particles between radii 
of 40 and 137 cm. The COT symmetry axis is the main 
axis of the cylindrical coordinate system used at CDF. 
Both tracking devices are immersed in a uniform axial 
1.4 T magnetic field, allowing precision measurement of 
the momenta of charged particles in the radial direction 
Pt- Planar drift chambers located outside of calorime- 
ters and additional steel absorbers are used to identify 
muons. Muons and charged hadrons are reconstructed in 
the central pseudorapidity region < 1. 

The three-level trigger system of CDF is used to select 
events of interest by requiring two oppositely charged 
particle tracks, each with pt > 1.5GeV/c and matched 
to hits in the muon detector. About two million J/tp —> 
/i + fi~ signal candidates were selected by requiring a re- 
constructed mass within 80.0 MeV/c 2 of the J/tp mass Q. 
A P° (B d ) meson candidate is reconstructed by associ- 
ating a J/tp candidate with a pair of tracks, each with 
Pt > 0.4GeV/c, consistent with a <p — ► K + K~ (K*° — > 
K + tt~) decay. The K + K~ (K + tt~) mass is required to 
be within 6.5MeV/c 2 (50.0MeV/c 2 ) of the 4> (K*°) pole 
mass Q. When both K+n- and tt + K particle assign- 
ments are kinematically viable in the K*° reconstruction 
(no particle identification is used), the one giving the 
mass closest to the pole K*° mass is chosen. This reduces 
the contribution from candidates with swapped particle 
assignment down to about 10% of the total signal. Back- 
ground is suppressed by requiring Pr{B a s d ) > 6.0GeV/c, 
p T {4>) > 2.0GeV/c and p T (K*°) > 2.6GeV/c. The more 
restrictive cut on px(K* ) is determined by an optimiza- 
tion procedure, and is a consequence of the larger com- 
binatorial background underneath the K*° peak. 

We fit the B candidates subject to the constraint that 
the four tracks originate from a common point. In order 
to improve the mass resolution, the mass is con- 

strained to the J/tp mass. To ensure only well measured 
vertices, a set of track and vertex quality requirements is 
applied 8j. In particular, all four tracks are required to 
have measurements in at least 3 axial layers of the silicon 
detector. The proper decay time, t, is determined from 
the radial distance It from the beam axis to the B me- 
son decay vertex, signed relative to the direction of pr 
of the B candidate: ct = c(It ■ Pt)Mb/Pt- The position 
of the beam axis is determined using data taken with an 
inclusive jet trigger. The radial profile of the beam is 
approximately Gaussian with an RMS of about 30 /an. 

In the B d system, the mass eigenstates are not CP 
eigenstates, and the observed decays are flavor-specific 
with the charge of the K meson identifying whether the 
decay is that of a B d or a B d . Summing over initially pro- 
duced B d and B d yields a differential decay rate which 
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is insensitive to first order to a lifetime difference 0: 
oc[| j 4 | 2 -/ 1 (a)+ L4 y | 2 -/ 2 (p) 



d 4 V(p,t) 



dpdt 



\A ± \ 2 -f 3 (p)±lm(AlA ± )-f 4 (p) 
-Re(^^n) • f 5 (p)±Im(A* A x ) ■ f 6 (p)]< 



-r d t 



Here the upper (lower) sign is used for K + ir (K ir + ) 
in the final state and Td = \/t b o. The functions fi(p) 
depend on the transversity variables p = {cos 6, tp, cos ?/>}, 
all of which are defined in Ref. |9(. 

In the B® system, the SM expectation is that CP vi- 
olation due to mixing is small, and the mass eigenstates 
are nearly CP eigenstates. Ignoring CP violation due to 
mixing, and summing the distributions for initially pro- 
duced and B s mesons, the interference terms between 
the CP-even and CP-odd amplitudes cancel, leaving: 



d 4 V(p,t) 
dpdt 
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/ 3 (p)+Re(^A||)-/ 5 (p) e 



-v L t 



In the P° analysis, the observed final states have definite 
CP which depends on L. The CP-even angular decay 
terms (A 0j ||) evolve in time as e~ TLt , while the CP-odd 
angular decay terms (A±) evolve as e~ Tllt . 

For both P° an d B° s decays, the amplitudes are nor- 
malized so that |A | 2 + \A\\ | 2 + \A±_\ 2 — 1, and an unob- 
servable overall phase is removed by setting arg(Ao) = 0. 
The decay amplitudes are assumed to be CP conserving. 

An unbinned likelihood fit is performed to mass, ct, 
and p in order to extract the decay amplitudes Aq im 
and decay widths T (Th and in the case of the B° s ). 
Inclusion of the mass information in the fit is crucial for 
separation of the signal from the background. The mass 
distribution is modeled with a Gaussian for the signal 
peak and a linear shape for the background. The mass- 
measurement uncertainty is incorporated for each can- 
didate. The probability density function for ct includes 
positive exponentials for the signal, a ^-function for the 
prompt background (which is about 85% of the total) 
and a set of exponentials for positive and negative decay 
lengths which describe a short-lived background compo- 
nent due to mis-measured vertices and a long-lived con- 
tribution due to incorrectly reconstructed heavy-flavor 
decays. Each contribution is convoluted with a Gaussian 
resolution function, the width of which is proportional to 
the uncertainty of the candidate's ct measurement. To 
allow for a systematic underestimate of the uncertainties, 
the mass and the ct uncertainties are multiplied by scale 
factors determined in the fit. The p distribution of the 
signal is parameterized in accordance with the equations 
above. The background distributions in ct and p are as- 
sumed to be uncorrelated. The latter is described by a 
shape similar to that of the signal, but with an indepen- 
dent set of amplitudes. The relationship of mass, ct, and 
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FIG. 1: Mass distribution with the fit projection overlaid: 
(a) B° 3 ^J/^t, (b) B° d -» J/i> K*° . 



p of the P° candidates with Kit mis-assignment to those 
of correctly reconstructed candidates is established via 
Monte Carlo simulation. 

Distributions in p are distorted by the detector accep- 
tance, the trigger efficiency and, most importantly, the 
kinematic selection criteria. We use the method devel- 
oped for the CDF Run I measurement of transversity 
amplitudes [ToL ITlj to account for this distortion. With 
as little as six constants extracted from Monte Carlo de- 
cays generated uniformly in p, this method allows one to 
avoid the need for explicit parameterization of the dis- 
tortion in the likelihood. All aspects of the fitting are 
extensively verified using Monte Carlo simulations. 

Data and fit projections in mass and ct for the B® and 
P" are shown in Figures Q]and Fits in the transversity 
sub-space are illustrated by Figure 

The single largest source of systematic uncertainty in 
the measurement of the transversity amplitudes of P° 
and B® is the choice of parameterization of the back- 
ground distribution in p. The B° s transversity amplitudes 
receive a small contribution to their systematic uncer- 
tainty from a 3.5% contamination from Two other 
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FIG. 2: ct distribution with the fit projection for the signal 
and background (bkg.) overlaid: (a) — > J/t/j<fr , (b) B® — > 
J/^K r0 . 



sources contribute to the uncertainty in £>° amplitudes: 
the way candidates with incorrect Kit assignment are 
handled and the potential contribution of p + p~ K + n~ fi- 
nal states that are not due to — > J/tpK*°(892) decays, 
which is estimated to be less than 4% of the total signal. 
The systematic uncertainty in the lifetimes receives con- 
tributions from the choice of the ct parameterization and 
from the SVX alignment. Slightly larger contributions 
come from the choice of the background parameterization 
in p and, in the case of B®, from B% contamination. For 
AT S /T S the only two sources of systematic uncertainty 
are from the choice of the background p parameterization 
and from B% contamination. Other potential sources of 
systematic uncertainty, including those from the method 
of correcting for distortion of the signal distribution in p 
and potential contribution of B° s — > J / 'tp f° (980) decays, 
were found to be negligible. The precision of all of the 
results of this analysis is statistically limited. 

Results from the time-dependent angular analysis are 
given in Table [I] The B% decay amplitudes and phases 
are of comparable precision and in agreement with re- 



TABLE I: Summary of the results of the time-dependent an- 
gular analysis of B° — > J/ip<j> and B° — > J/tpK*° decays. A 
measurement of arg(Aj_) is not possible for the B® decays be- 
cause the final state particles do not distinguish the decay of 
a from that of a B s . For the B® decays, any pair of quan- 
tities describing signal lifetimes (t b o = 1/T S , th, t~l, AT 3 /r s 
and AT S ) may be used as free parameters in a fit; separate 
fits using different pairs are performed to obtain directly the 
results and asymmetric (statistical) uncertainties for each of 
the quantities. 



B' 



N stg 

A 

I All 
\A±\ 

ar g(^||) 
arg(A ± ) 



203 ± 15 

0.784 ± 0.039 ± 0.007 
0.510 ± 0.082 ±0.013 
0.354 ±0.098 ±0.003 
1.94 ± 0.36 ± 0.03 



CTgO 


(419 til ± 6) pm 


CT B 


CTL 


(316 ttl ± 6) pm 


CT H 


(622 till ± 9) pm 


Ar s /r s 


(65 tf 3 ± 1) % 


Ar s 


(0.47 t° ;ll ± 0.01) ps 



1155 ±39 

0.750 ±0.017 ±0.012 
0.473 ± 0.034 ± 0.006 
0.464 ± 0.035 ± 0.007 
2.86 ±0.22 ±0.07 
0.15 ±0.15 ±0.04 

(462 ± 15 ± 6) pm 



suits from BaBar and Belle and the lifetime is 
in agreement with the world average value 0. Previous 
results [7l]] for the B® decay amplitudes, obtained from 
a time- integrated analysis, are in agreement with the re- 
sults obtained in this analysis. Within uncertainties, the 
amplitudes for the B° and B® decays are in agreement, 
as is expected in the limit of SU (3) flavor symmetry [j| . 
Explicitly requiring exact SU(3) symmetry by setting the 
B° decay amplitudes to be equal to those of the S° , gives 
a consistent result for Ar s /T s within uncertainties. 

It is predicted [3j that the B^ and B® total decay 
widths should be equal to within 1%. This expecta- 
tion can be used as a constraint in the B® fit by re- 
quiring l/r s = T B = 2T H T L /(T H + T L ) = T B EE 1/T d , 

with ct b o = 460.8 ± 4.2 pm, the known value for the J5° 
lifetime with an additional 1% uncertainty added in 
quadrature. By applying this constraint in the fit, we find 
i and Ar s = (0.46 t° ,ll ± 0.01) 
1 Although the uncertainties are still sizable, the fits 



AT./r, = (71 +_lt ± I 
ps 1 

with and without the T s — constraint favor a non-zero 
value for Ar s /r s . 

Monte Carlo methods are employed to estimate the 
probability for an experiment with similar statistical sen- 
sitivity to yield Ar s /r s as large as is observed in this 
analysis. For the SM expectation, Ar s /r s = 12%, one 
experiment in 84 (204) would give a result larger than 
that obtained from the unconstrained (constrained) fit. 
If no lifetime difference were expected, Ar s /T s = 0, one 
experiment in 315 (718) would give a result larger than 
that obtained from the unconstrained (constrained) fit. 



7 




A lifetime difference should result in an increase of the 
fraction of CP-odd B® — > decays obtained in a 
time-integrated angular fit as a function of a cut on ct. 
The time- integrated CP-odd fraction extracted from fits 
for four values of the ct cut are shown in the second col- 
umn of Table [n] Using ctx and cth obtained from the 
time-dependent fit and the observation that the time in- 
tegrated CP-odd fraction is 20% for ct > 0, one obtains 
expected fractions for the other ct cuts, shown in the 
third column in Table [H] which are in agreement with 
the observations. The P° — > J/ipK*° decays provide 
an important cross check of the results obtained for the 
B® — > J/tp <fi decays. Fits for the time-integrated fraction 
of parity-odd P° decays show, as expected, that the frac- 
tion remains unchanged with respect to cuts on ct (last 
column of Table ITl)l . In addition, a fit of the B^ data can 
be performed allowing two lifetime components. The re- 
sults are consistent with no lifetime difference for the full 
sample (AT/T = (15 ± 12)%), as well as for independent 
subsamples having a statistical sensitivity similar to the 
B® decay sample. This result is not a measurement of a 
lifetime difference in the B® system, but rather a cross 
check of the analysis technique. 

In conclusion, we have performed the first time- 
dependent angular analysis of B® — > J/tp <\> decays and 
have performed a similar measurement with B® — > 



TABLE II: Time-integrated CP-odd P° and parity-odd B\ 
fractions (in %) vs. a cut on the decay length, ct. 



ct cut 


B" fitted 


B" expected 


B^ fitted 


> fj,m 


20 ±9 


20 (reference) 


22 ±4 


> 150 


24 ± 10 


24 


23 ±4 


> 300 


30 ± 13 


29 


23 ±4 


> 450 /im 


39 ± 12 


34 


24 ±5 



J/ipK*° decays. The measured P° polarization am- 
plitudes are of comparable precision to, and in agree- 
ment with, previously published results. The measured 
B® polarization amplitudes are the most precise avail- 
able. Analysis of the P° — » J/ip<j> decays indicates a 
non-zero lifetime difference between the heavy and light 
mass eigenstates of the P° system. The result obtained, 
AF s /r s = (65 ± 1)%, has a central value larger than 
the SM expectation of (12 ± 6)%. 
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